Possible nodeless superconductivity in the noncentrosymmetric superconductor Mg^-sIr^E^ 
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We measured the resistivity, diamagnetization, and low-temperature specific heat of the newly discovered 
noncentrosymmetric superconductor Mg^-jIrigBig. The temperature dependence of specific heat is consistent 
with the model of an isotropic s-wave gap with value Ao ~ 0.94 meV for the sample T c = 5.7 K, and the ratio 
Ao/k B T c » 1.91 indicates a slightly moderate coupling for the superconductivity. The correlations among the 
normal state Sommerfeld constant y„, the slope -dpaH c i{T)ldT near T c , and the condensation energy E c are 
all consistent with the slightly moderate coupling picture. Based on the data of phonon contribution, T c and 
the McMillan formula, we obtained an electron-phonon coupling strength A e - P h ~ 0.66, which suggests that the 
superconductivity here is induced by the electron-phonon coupling. 

PACS numbers: 74.20.Mn,74.20.Rp, 74.25.Bt, 74.70.Dd 



I. INTRODUCTION 



II. SAMPLE PREPARATION AND CHARACTERIZATION 
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The study on superconductivity in noncentrosymmetric ma- 
terials has attracted growing efforts in recent year s 1 i 2 i 3 ' 4 ' 5 . For 
most superconductors, the atomic lattice has a centrosymme- 
try, therefore the system is inversion symmetric. The orbital 
part of the superconducting order parameter has a subgroup 
which is confined by the general group of the atomic lattice. 
Due to the Pauli's exclusion rule and the parity conservation, 
the Cooper pair with orbital even parity should have anti- 
parallel spin state, namely spin singlet, while those having or- 
bital odd parity should have parallel spin state, i.e., spin triplet. 
If a system lacks the centrosymmetry, the enhanced spin- 
orbital coupling may promote the formation of pairing with 
high angular momentum, such as the triplet pairing. Mean- 
while the noncentrosymmetric structure allows for the exis- 
tence of a mixture of singlet and triplet pairing. Theoretical 
features are anticipated in the noncentrosymmetric system^. A 
nodal gap structure has been observed in Li2Pt3B showing the 
possibility of triplet pairing, while due to weaker spin-orbital 
coupling 7 '^, the nodal gap has not been observed in a mate- 
rial Li2Pd3B with similar structure. It is thus highly desired to 
investigate the paring symmetry in more materials with non- 
centrosymmetric structure. 

The newly discovered superconductor MgioIri9Bi6 (here- 
after abbreviated as MglrB ) (Ref. 9) with superconducting 
transition temperature T c ~ 5 K is one of the rare materi- 
als which have the noncentrosymmetry. This material has 
a space group of /-43m with large and complex unit cells 
of about 45 atoms. To some extent it resembles the system 
Li2(Pt,Pd)3B since they have alkaline metals (Li, Mg), heavy 
transition elements (Pd, Pt, Ir) and the light element boron. 
Theoretically it was shown that the major quasiparticle den- 
sity of states (DOS) derives from the d orbital of the heavy 
transition elements. In this paper we present a detailed inves- 
tigation and analysis on the superconducting properties, such 
as the energy gap, pairing symmetry, electron-phonon cou- 
pling strength, condensation energy, etc., in MglrB. Our re- 
sults suggest that the superconductivity in this system is of 
BCS type with an s-wave gap symmetry and a slightly mod- 
erate electron-phonon coupling strength. 



The samples were prepared in two steps starting from pure 
elements of Mg (98.5%), Ir (99.95%) and, B (99.999%) us- 
ing a standard method of solid state reaction. Appropriate 
mixtures of these starting materials were pressed into pellets, 
wrapped in Ta foil, and sealed in a quartz tube with an atmo- 
sphere of 95% At/ 5% H2. The materials were then heated 
at 600 °C and 900 °C for 40 min and 80 min, respectively. 
After cooling down to room temperature, the samples were 
reground and then they were pressed into pellets and sealed 
in a quartz tube with the same atmosphere as used in the first 
step. Some of the time they were mixed with another certain 
amount of Mg up to 20%. In this process the sample was 
heated up to 900 °C directly and maintained for 80 min. Usu- 
ally the superconducting transition becomes sharper after the 
second process, but sometimes the T c may become slightly 
lower than the first time. The synthesizing process here is 
similar to the previous work reported by the Princeton group 9 
but still with some differences. For example, we used Mg 
powder instead of flakes to make the mixture more homoge- 
neous. In addition, the pressure in the sealed quartz tube may 
rise to nearly 4 atm at 900 °C, which may considerably re- 
duce the volatilization of Mg during the synthesis. This is 
probably the reason that the sample with the transition tem- 
perature as high as T c — 5.7 K was made without adding extra 
Mg in the second step with the starting ratio in the first step as 
Mg:Ir:B=12:19:16. 

The resistivity and the ac susceptibility were measured 
based on an Oxford cryogenic system (Maglab-Exa-12). The 
specific heat was measured on the Quantum Design instru- 
ment physical property measurement system (PPMS) with 
temperature down to 1.8 K and the PPMS based dilution re- 
frigerator (DR) down to 150 mK. The temperatures of both 
systems have been well calibrated showing consistency with 
an error below 2% in the temperature range from 1.8 K to 10 
K. 

As shown in Fig.l, the x-ray diffraction (XRD) pattern 
taken on one sample with T c — 5.7 K shows a single phase 
with very small amount of impurity which is comparable to 
that reported previously 9 . After the first round of synthesiz- 
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FIG. 1 : (Color online) The x-ray diffraction pattern measured for the 
sample T c = 5.7 K. The peaks from the secondary impurity phase are 
marked by the arrows. It is clear that the main diffraction peaks are 
from the phase MglrB. 
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FIG. 2: (Color online) Temperature dependence of resistivity (top) 
and magnetic susceptibility and^f') (bottom) under different dc 
magnetic fields. It is clear that the dc magnetic field makes the tran- 
sition shift parallel to low temperatures, manifesting a field-induced 
pair-breaking effect. The inset in the top panel shows the resistive 
transition in a wide temperature regime at fi H = T. 



ing, the superconducting transition inspected by the ac suscep- 
tibility occurs at about 5 K with a relatively wide transition. 
However, after the second step, the transition moves to about 
5.7 K with a sharper transition width if we did not add extra 
Mg in the second step. The Ta foil seems unreacted with the 
materials in both steps of the fabrication if the Mg content is 
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FIG. 3: (Color online) Raw data of specific heat plotted as C/T vs 
T 2 . All filled symbols represent the data taken with PPMS at various 
magnetic fields. The open symbols show the data taken with the DR. 
The thick solid line represents the normal state specific heat which 
contains both the phonon and the electronic contributions. 



below 12 in the stoichiometric formula. 

In the top frame of Fig. 2 we show the temperature de- 
pendence of resistivity under different magnetic fields. One 
can see that the transition width determined from resistive 
measurements (1% - 99%p„ ) is only about 0.2 K. This 
is consistent with the rather sharp magnetic transition as 
revealed by the ac susceptibility data shown in the bot- 
tom view of Fig. 2. By applying a magnetic field the 
transition shifts to lower temperatures quickly with a slope 
-d/JoH C 2(T) I dT\r c ~ 0.63 T / K for the present sample with T c 
= 5.7 K. Using the Werthamer-Helfand-Hohenberg relation 10 
HqH C 2 - -0.69dfXoH C 2(T) /dT\r c T c , we get the upper critical 
field fj.oH C 2 = 2.48 T. The value found here is comparable to 
that reported in the earlier paper 9 (T c * 5 K). It is interesting 
to note that the value of the slope -dfioH C 2(T)/dT correlates 
strongly with T c . For a sample with T c = 3.7 K, we found that 
-dfj.()H C 2{T) I dT\r c ~ 0.2 T / K. This suggests that there is a 
tunability of superconducting properties in this system which 
will be addressed in Sec. IV. 



m. SPECIFIC HEAT AND FIT TO THE BCS MODEL 

Shown in Fig. 3 are the raw data of the specific heat up 
to a magnetic field of 1 T. The open symbols represent the 
data taken with the DR, while all filled symbols show the data 
taken with the PPMS. Both sets of data coincide very well for 
different fields. With increasing the magnetic field the spe- 
cific heat jump due to the superconducting transition moves 
quickly to lower temperatures leaving a background which is 
consistent with that above T c at zero field. This provides a 
reliable way to extract the normal state specific heat as shown 
by the thick solid line since the normal state can be described 
by C/T — y n +/3T 2 , where the first and the second term corre- 
spond to the normal state electronic and phonon contribution, 
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DOS at Ef given by the local-density approximation band- 
structure calculation 11 is about N(Ef) = 5.51/eV spin. By 
assuming that there is an electron-phonon coupling constant 
A e -ph in the system, one has 



2tt 2 , 

Jn = —N(E F )k 2 B (l + Ae-pk). 



(1) 



Using the band-structure value of N(Ff), we have y n - 
25.98 (1 + /l e _ p /,)(mJ/mol K 2 ). Taking the experimental value 
y„ = 43.9 mJ/molK 2 , we obtain A e -ph ~ 0.68, indicating a 
slightly moderate electron-phonon coupling. 

In the raw data shown in Fig. 3, one can see an upturn of 
y = C/T in the very low-temperature region. This upturn is 
known as the Schottky anomaly, induced by lifting the degen- 
eracy of the states of the paramagnetic spins. We tried a two 
level (S=l/2) model to fit the low-temperature data but found 
a poor fitting together with an extremely large Lande factor g 
in the Zeeman energy gfj.BH e fj, where /ub is the Bohr magne- 
ton, H e ff - tJh 2 + Hq is the effective magnetic field which 
evolves into = Hq, the crystal field at zero external field. 
In MglrB the most possible paramagnetic centers may be from 
Ir 4+ (S=5/2) or Ir 3+ (S=2). The system energy due to Zeeman 
splitting in a magnetic field is 12 

E S ch = 2 E i e*p(-Ei/k B T)/ ^ exp(-Ei/k B T), (2) 



FIG. 4: (Color online) Temperature dependence of y e - y„ for (a) 
magnetic fields up to 1 T and (b) at zero field before (filled circles) 
and after removing (open and filled squares) the Schotkky anomaly. 
The dark solid, red dashed, and dotted green lines in (b) are the the- 
oretical curves calculated based on the BCS model (see text) with a 
gap of s wave, d wave, and p wave, respectively. 



respectively. From the data it is found that /3 — 5.4 mJ/ mol 
K 4 and y„ = 57.7 mJ/mol K 2 . By extrapolating the data in 
the superconducting state at zero field down to K one finds, 
however, a residual value yo * 13.8 mJ/ mol K 2 indicating a 
contribution either by a nonsuperconducting fraction in vol- 
ume of about 24%, or by the DOS induced by the impurity 
scattering for a nodal gap which will be discussed later. The 
residual yo may be induced by the nonsuperconducting region, 
but it is still difficult to be regarded as due to an impurity phase 
with completely different structure as MglrB since the XRD 
data shown in Fig. 1 is quite clean. We thus suggest that the 
superconductivity depends sensitively on the relative compo- 
sitions among the three elements and some regions without su- 
perconductivity have the chemical composition and even the 
structure close to the superconducting phase. In any case, it is 
safe to conclude that the normal state Sommerfeld constant for 
the present sample is close to or slightly above 43.9 mJ/mol 
K 2 . 

Next we can have an estimation on the electron-phonon 
coupling strength based on y n . In MglrB, the electronic con- 
duction is dominated by the 5d band electrons of Ir atoms. The 



where E t = Mjgfi B H ef f and Mj = -S, -S+l, S-l, S. 
The specific heat due to the Schottky effect is thus Cs c h = 
(n/kB)dEsch/dT, where n represents the concentration of the 
paramagnetic centers. For 5=5/2 and S — 2 the calculated 
results are very close to each other, therefore we show only 
the fit with S = 5/2 corresponding to Ir 4+ (six levels). This 
method allows us to deal with the data at zero and finite fields 
simultaneously. It is known that the Schottky term should be 
zero at T — K. In the superconducting state, the total spe- 
cific heat can be written as C tot = C nons + C e + C p h + Cs c h 
with Cnons = yoT as the contribution of the nonsuperconduct- 
ing regions, C e is the electronic part. In the zero temperature 
limit only the contribution of the nonsuperconducting part is 
left. Applying a magnetic field gives rise to a finite value Ay e 
to C e = y e T due to the presence of vortices. Practically, in or- 
der to fit the Schottky term, we first remove the phonon con- 
tribution Cph = /3T 3 , then vertically move the experimental 
data downward with a magnitude yo = 13.8 mJ/ mol K 2 and 
a field-induced vortex term Ay e (H). In Fig. 4(a) we present 
the specific heat coefficient y e - y„ measured at the magnetic 
fields up to 1 T. An upturn due to the Schottky effect is visible 
in the low temperature region for all fields. In Fig. 4(b) we 
show the data at fi^H — T before [(red) filled circles] and 
after (dark open squares) removing the Schottky anomaly cal- 
culated based on the model of six levels with h^Hq = 0.25 T, 
n = 3.74 mJ/mol K and g = 2. Note that we used the value 
43.9 mJ/ mol K 2 as the normal state Sommerfeld constant y„. 
One can see that the low-temperature part after removing the 
Schottky anomaly is flattened out below about 0.8 K when the 
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field is zero. Furthermore, it can also be justified by the re- 
quirement of entropy conservation. Since the Schottky term 
gives only a very small contribution in the high-temperature 
region (above 1 .5 K here), if y e had a power law as required 
by a nodal gap, instead of a flat temperature dependence for 
an i-wave gap, the entropy would be clearly not conserved 
yielding a large negative entropy. This is of course unreason- 
able. In Fig. 4(b) we present together the theoretical curves 
for y e - y n calculated using the weak coupling BCS formula 
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where f = ye 2 + A 2 (T, 0). In obtaining the theoretical fit 
we take the implicit relation Aq(T) derived from the weak 
coupling BCS theory for superconductors with different pair- 
ing symmetries: A(T,9) = Aq(T) for s wave, A(T,ff) - 
Aq(T)cos28 for d wave, and A(T, 0) = Aq(T)cos6 for p wave, 
respectively. The theoretical curve of s wave fits the experi- 
mental data very well leading to an isotropic gap value Ao = 
0.94 meV and T c = 5.7 K. The ratio A /k B T c = 1.91 obtained 
here is quite close to the prediction for the weak coupling 
limit (Ao I 'IcbT c = 1.76), indicating a slightly moderate cou- 
pling strength. This is self-consistent with the conclusion de- 
rived from the estimation on y„. In addition, the specific heat 
anomaly at T c is AC e /y n T\r c ~ 1 .64 being close to the theoret- 
ical value 1 .43 predicted for the case of weak coupling, again 
showing a slightly moderate coupling. The inset in Fig. 4(b) 
shows a field induced part Ay e . In an s-wave superconductor 
Ay e is mainly contributed by the vortex cores and a linear rela- 
tion Ay e cc Hy„/H C 2(0) is anticipated 13 in the low field region 
with AZ/Ep <k T <k T c . This linear relation is well demon- 
strated by the data below 0.7 T, indicating another evidence of 
i-wave pairing symmetry. This is in sharp contrast with the re- 
sults in cuprates where a Ay e oc yfH relation is expected due 
to the Doppler shift of the nodal quasiparticle spectrum in a 
c/-wave superconductoriii^. The theoretical curves calculated 
with af-wave and p-wave models cannot describe the data in 
both the low-temperature and high-temperature regime s 13 ' 16 . 

In the following we try to estimate the superconducting 
condensation energy E c for the sample with T c = 5.7 K. In 
calculating E c we get the entropy difference between the nor- 

mal state and the superconducting state by S „ - S s — J Q (y„ — 

y e )dT', then E c is calculated through E c = Jj K (S„ - S s )dT' . 
The data of S„ and S s as well as the difference between them 
are shown in insets (a) and (b) of Fig. 5, respectively. The 
main frame of Fig. 5 shows the temperature dependence of 
the condensation energy E c which is about 369 mJ/mol at T 
= OK. This value can actually be assessed by the following 
equation 



E c = aN(E F )A 2 /2 = a ^^y„A 2 . 
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FIG. 5: (Color online) The main frame shows the superconduct- 
ing condensation energy calculated from the specific heat. Inset (a) 
shows the entropy in the normal and superconducting state. Plotted 
in inset (b) is the difference of the entropy between the normal and 
superconducting state. 
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For a BCS s-wave superconductor, a = 1, taking y„ = 43.9 
mJ/mol K 2 and Ao = 0.94 meV, we found a value of E c as 
367 mJ/mol which is remarkably close to the experimental 
value 369 mJ/mol. This also validates the values of Ao and y„ 
determined in our experiment. 

Now we discuss the possibility of impurity scattering which 
may suppress the spin-triplet component of the superconduc- 
tivity. Theoretically it is known that the nodal gap structure is 
very sensitive to the impurities. If the spin-singlet and triplet 
components are mixed, the latter might be suppressed by the 
impurity scattering and the system would behave like a BCS 
superconductor. Although we cannot exclude this possibility 
at this moment, it is safe to conclude that the spin-triplet com- 
ponent (if it exists) should be a very small part of the total 
condensate. This can be justified by the following two ar- 
guments. First, for the sample with T c = 5.7 K, we found 
that yo/y„(T > T c ) = 24%. Assuming that the total yo were 
induced by the impurity scattering of this spin-triplet compo- 
nent, the ratio of the DOS corresponding to the normal state 
of this part would be less than 24%. Actually one would 
not expect that the total value of yo here could be ascribed 
completely to the suppressed spin-triplet component since the 
XRD data does show the existence of some secondary (non- 
superconducting) phase which will certainly contribute a cer- 
tain value of DOS. Second the self-consistency among the de- 
rived values of the normal state Sommerfeld constant, the ap- 
propriate ratio of Au/IcbT c and the electron-phonon coupling 
strength (see next section) together with the nice fit of the spe- 
cific heat data to the BCS model all can push the spin-triplet 
component to a very low limit. Nevertheless, as a future work 
it is worthwhile to have an inspection on this issue on a sample 
with complete purity. 
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FIG. 6: (Color online) Temperature dependence of resistivity (top) 
and magnetic susceptibility and^') (bottom) under different dc 
magnetic fields for the sample with T c = 3.7 K. It is clear that the dc 
magnetic field makes the transition shift parallel to low temperatures 
with a much faster speed compared with the sample with T c = 5.7 K. 
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FIG. 7: (Color online) The upper critical field for both samples with 
different transition temperatures. 



IV. ELECTRON-PHONON COUPLING STRENGTH AND 
THE TUNABILITY OF SUPERCONDUCTIVITY IN MGIRB 

Now we get to the electron-phonon coupling in MglrB. 
From the normal state value we have derived the phonon con- 
tribution C ph = J3T 3 with j6 « 5.4 mJ/mol K 4 . Using the re- 



lation @ D = (\2n' i k B N A ZI5l3y l \ where N A = 6.02 xlfy j is 
the Avogadro constant, Z = 45 is the number of atoms in one 
unit cell, we get the Debye temperature D (MgIrB) = 253 K 
for the sample with T c — 5.7 K. If the conduction electrons 
are weakly coupled with the phonon, one can use the McMil- 



lan equation to evaluate the electron-phonon coupling strength 

X e - ph viaii 



®D 

1.45 



exp[- 



1.04(1 + A e - ph ) 



H*{\+0.62A e - ph ) 



U-ph 



(5) 



where yt is the Coulomb pseudopotential taking 0.11. Us- 
ing ©o(MglrB) = 253 K and T c = 5.7 K, we obtain A e - ph 
= 0.66, indicating a slightly moderate coupling strength. The 
value determined here is quite close to A e - P h = 0.68 which was 
determined previously based on the comparison between the 
normal state Sommerfeld coefficient y„ and the theoretically 
estimated DOS at £>. All these detailed and self-consistent 
analysis indicate an electron-phonon coupling mechanism for 
superconductivity in MglrB. 

At this moment, it is still unclear that the conduction elec- 
trons are strongly coupled to which phonon branch leading to 
the superconductivity. During the preparation of samples, it 
was found that the superconducting transition temperature is 
changeable in a wide region (from 2 K to 5.7 K) in MglrB with 
a rather sharp transition width. This is actually a rare case in 
the alloy superconductors in which the off-stoichiometric el- 
ements normally play as the scatterers and break the Cooper 
pairs and finally make the transition broad. In MglrB a slight 
off-stoichiometry may not change the structure, but promote 
the mutual-substitution leading to a different electron-phonon 
coupling strength. In Fig. 6 we show the temperature de- 
pendence of the resistivity and the ac susceptibility of another 
sample (T c = 3.7 K) which was made also in two steps but 
with 20% more Mg added in the second step. One can see 
that the transition width at zero field is still quite narrow and 
the x-ray diffraction pattern taken on this sample is reason- 
ably clean (not shown here). An estimation on the electron- 
phonon coupling strength A e - P h based on the McMillan equa- 
tion on this sample tells that ~ 0.55. In Fig. 7 we 
present the upper critical field H C 2(T) taken with the crite- 
rion 10%p„ for the two selected samples with T c = 5.7 K 
and T c = 3.7 K. One can see that the slope -d^ioH C 2/dT\ Tc 
is very different between these two samples. Actually in the 
previous work of the Princeton group 9 , -dii( l H C 2ldT\j c can 
be as high as 1 T/K. In a dirty type-II superconductor, it was 
predicted 18 that -dfioH C 2/dT\r c °c p n Jnr\ with r\ in connection 
with the electron-phonon coupling strength. It is thus reason- 
able to ascribe the variation of d/ioH C 2/dT\T c and T c to differ- 
ent electron-phonon coupling strengths and different density 
of states at £>. All these quantities may be optimized in fu- 
ture work and hopefully will lead to a higher superconducting 
transition temperature. The basic parameters and properties 
derived in this work provide a playground for the future study 
in this interesting system. 



V. CONCLUDING REMARKS 

In summary, analysis on the low-temperature data in MglrB 
finds an s-wave pairing symmetry with a gap in the weak 
coupling to slightly moderate coupling regime and the su- 
perconductivity is induced by the electron-phonon coupling. 
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The spin-triplet component, if it exists, should remain as a 
small part of the total condensate. Tuning the superconducting 
transition temperature seems possible through changing the 
electron-phonon coupling strength and the density of states at 
the Fermi level by varying the relative compositions among 
the three elements and probably also by mutual substitution. 
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